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We report the first lattice QCD calculation of the complex kaon decay amplitude Ao with physical 
kinematics, usinga32^x64 lattice volume and a single lattice spacing o, with 1/a = 1.3784(68) GeV. 
We find Re(Ao) = 4.66(1.00)(1.26) x 10“^^ GeV and Im(Ao) = —1.90(1.23)(1.08) x 10“^^ GeV, where 
the first error is statistical and the second systematic. The first value is in approximate agreement 
with the experimental result: Re(Ao) = 3.3201(18) x 10“^ GeV while the second can be used to 
compute the direct GP-violating ratio Re(£'/£) = 1.38(5.15)(4.59) x 10“'^, which is 2.1 (t below the 
experimental value 16.6(2.3) x 10“^. The real part of Ao is GP conserving and serves as a test of 
our method while the result for Re(£'/£) provides a new test of the standard model theory of GP 
violation, one which can be made more accurate with increasing computer capability. 

PACS numbers: 12.38.Gc ll.30.Er 12.15.Hh 13.20.Eb 


The violation of CP symmetry was discovered as a sub¬ 
percent admixture of the CP-even combination of K® 
and K® mesons in a nominally CP-odd decay eigen¬ 
state [ij. In the standard model this mixing is caused 
by a single CP-violating phase which can be introduced 
if there are three generations of quarks in nature 0 . This 
CP-violating mixing is the indirect effect of virtual top 
quarks. It is described by the parameter e whose mea¬ 
sured magnitude is 2.228(0.011) x 10“®, a value success¬ 
fully related by the standard model to the CP-violating 
phase measured in the decay of bottom mesons. 

Much more difficult to measure and to compute the¬ 
oretically is the direct violation of CP in K decay, de¬ 
scribed by the parameter e' and resulting from a CP- 
violating difference between the phases of the decay am¬ 
plitudes Aq and A 2 , which describe kaon decay into 
a two-pion state with isospin 1 = 0 and 2, respec¬ 
tively. This direct CP violation is 3 orders of magnitude 
smaller than that caused by mixing, with Re(e'/£) = 
1.66(0.23) X 10-® Because of its small size this di¬ 
rect violation of CP is especially sensitive to phenomena 
beyond the standard model, phenomena that are believed 
to be required to explain the current excess of matter over 
antimatter in the Universe. 


in QCD and electroweak perturbation theory 


Hw = ^V*,yud'^[zi{T) + ( 1 ) 

^ i—1 

Here Gp = 1.166 x 10“®/(GeV)^, Vq'q is the Cabibbo- 
Kobayashi-Maskawa matrix element connecting the 
quarks q' and q and r = — V)* V)d/ 14 *gt 4 d- The ten opera¬ 
tors Qi are combinations of seven independent four-quark 
operators renormalized at the scale y.. The task that 
remains is to compute the matrix element of the ten Qi 
between an initial kaon and final tttt state with I = 0 or 2. 
While this has been an active area for theoretical work 
over the past thirty years, no reliable analytic method 
to compute these matrix elements has emerged Ifil-1^. 
However, this task is well suited to lattice QCD. 

Over the past five years, the calculation of the 1 = 2 
decay has become accessible to lattice methods [3, EBl 
and physical, continuum-limit results for A 2 are available 
with 10% errors [3]- However, calculating the 1 = 0 
amplitude Aq faces substantial new difficulties: i) the 
need to create an / = 0 two-pion state with energy well 
above threshold and ii) the statistical noise associated 
with the vacuum intermediate state. These difficulties 
have been overcome by methods we will now describe. 
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While standard model, direct CP violation involves 
massive W bosons and top quarks at an energy scale far 
above that accessible to lattice QCD, these high-energy 
interactions can be accurately captured by a low-energy 
effective Lagrangian with Wilson coefficients {yt and Zi 
below) which have been computed to next-leading-order 


COMPUTATIONAL METHOD 

The K —>■ TTTT matrix elements of the ten operators Qi 
are determined from the Euclidean Green’s functions 

CK,rritK,tQ, ^TTTT ) = {0\J^^{U^)Q,itQ)JK{tK)\0) (2) 
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in the limit of large time separations t,rTr ~ tg and tg — tx 
which projects onto the initial and final states of interest. 
The operators Jk and create the initial-state kaon 
and destroy the two final-state pions. Introducing a final 
state composed of two pions with nonzero relative mo¬ 
mentum poses special challenges. Using now standard 
methods 17|, the desired finite-volume two-pion state 
would have an energy well above that with two pions 
at rest and require a multiexponential fit to determine 
the decay matrix element. For the 1 = 2, two-pion state 
this problem can be addressed by imposing antiperiodic 
boundary conditions on the down quark Bin. 

However, for the 1 = 0 state we must impose isospin- 
symmetric boundary conditions to avoid mixing the I = 
0 and 2 states. This is possible through a major algo¬ 
rithmic advance: the introduction of G-parity boundary 
conditions Since each pion is odd under G par¬ 

ity, apart from the effects of their interaction, each pion 
must then carry a minimum momentum of tt/L for each 
direction (of length L) in which G parity is imposed. For 
our lattice volume, imposing G-parity boundary condi¬ 
tions in all three spatial directions results in the required 
/ = 0, TTTT energy ~ Mk- 

The G-parity transformation is described by the opera¬ 
tor G = Ge“'^“, a product of charge co njug ation (G) and 
a 180° isospin rotation about the y axis [^. When a lat¬ 
tice derivative connects quark fields across such a bound¬ 
ary the {u, d) doublet is joined to a G-parity transformed 
doublet (d, —u). This doubles the computational cost 
and requires substantial code modifications since explicit 
u and d degrees of freedom must be introduced. In addi¬ 
tion, the gauge helds must now obey charge-conjugation 
boundary conditions which demands new, special, gauge 
ensembles. Since quarks and antiquarks are mixed at the 
boundaries, new contractions must be included in which 
two quark or two antiquark fields are joined by a prop¬ 
agator. Finally, a consistent treatment of the strange 
quark s requires that we include an unphysical partner 
s' to form an isodoublet that obeys G-parity boundary 
conditions (2^ . When generating the 2-1-1 flavor gauge 
ensemble we must then take the square root of the deter¬ 
minant of the s — s' Dirac operator so that only a single 
strange quark flavor is included. 

The second critical difficulty is that the / = 0, two-pion 
state has the same quantum numbers as the vacuum, the 
state which thus dominates the large Gtt— tg limit needed 
to remove excited states. We must subtract this vac¬ 
uum contribution and deal with the exponentially falling 
signal-to-noise ratio that results, a subtraction carried 
out successfully in threshold calculations, with final-state 
pions approximately at rest. 23|, [2^ 25l |. 

We reduce the noise from this vacuum subtraction 
using two techniques. First, we use a split-pion op¬ 
erator B to destroy the two-pion state. Specifically, 
is the product of two quark-antiquark pairs, one 
pair at the time Gtt and the second at G7r+4. By separat¬ 


ing the pion operators we suppress the vacuum coupling 
that results when coincident pion operators immediately 
create and destroy a pion, reducing the vacuum noise 2 x. 
Second, we use all-to-all propagators to construct 

each pion interpolating operator from a quark-antiquark 
pair, fixed to Coulomb gauge, with a relative coordinate, 
hydrogen ground-state wave function of radius 2a and 
center-of-mass coordinate distributed over a time plane 
at t^r-TT or + 4. This choice increases the coupling 
to the two-pion state relative to the vacuum, giving a 
further 2x noise reduction [2^ . 

We use a 32^ x 64 volume, the IwasakU-DSDR gauge 
action [29| and Mobius [30|, domain wall fermions 
(DWF) [3]J with an extent of 12 in the fifth dimension. 
By using /3 = 1.75 and Mobius parameters 6 -|- c = 32/12 
and 6 — c = 1 we ensure that this ensemble is equivalent 
to our earlier dislocation-suppressing determinant ratio 
(DSDR) ensemble [32|, except that the latter has periodic 
boundary conditions and = 170 MeV. Input quark 
masses of m;(= = rud) = 0.0001 and nis = 0.045 

are used. (If a dimensioned quantity is given without 
units, lattice units are implied.) The inverse lattice spac¬ 
ing, residual quark mass, pion mass, and singlc-pion en¬ 
ergy are 1/a = 1.3784(68) GeV, = 0.001842(7), 
= 143.1(2.0) MeV and = 274.6(1.4) MeV. 

We analyzed 216 gauge configurations separated by 
four units of molecular dynamics time, starting at 300 
time units for equilibration. Seventy-five distinct dia¬ 
grams were computed, of four types as shown in Fig. [T] 
We compensated for this small number of configurations 
by performing 64 measurements on each conhguration, 
introducing the kaon and pion sources on each of the 
64 time planes. (The statistically more accurate, type 
1 and 2 diagrams were computed only on every eighth 
time plane.) The many propagator inversions needed 
on each configuration were accelerated usiM low-mode 
deflation with 900 Lanczos eigenvectors [33[ with the 
BAGEL fermion matrix package 3^. A complete set of 
measurements required 20 hours on an IBM Blue Gene/Q 
i-rack 33 , in balance with the 24 hours needed to gen¬ 
erate four time units of gauge field evolution on this same 
machine. 



types type4 


FIG. 1. Examples of the four types of diagram contributing to 
the AI = 1/2, K —y tttt decay. Lines labeled I or s represent 
light or strange quarks. Unlabeled lines are light quarks. 
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We must deal with two sorts of finite-volume effects. 
The first are errors falling exponentially with increasing 
lattice size which result from “squeezing” the physical 
states. Such errors are at the percent level if > 4. 
In our case, = 3.2 and errors ~ 7% may result [l5| . 
The second are effects falling as a power of L, similar 
to the discretization of the energy that we are exploit¬ 
ing. Here we apply the Lellouch-Luscher correction 0 
to remove the leading 1/L^ effect. This requires that our 
final TTTT state is an “s-wave” combination of the eight 
single-pion momenta (±1, ±1, ±l)7r/L. Ensuring this s- 
wave symmetry requires pion operators constructed to 
minimize the quark-level, cubic-symmetry violations in¬ 
troduced by G-parity boundary conditions. 

Essential to this calculation is the ability to define the 
seven independent, four-quark, lattice operators which 
correspond to those in the continuum Eq. (IT|). This is 
accomplished by using DWF whose accurate chiral sym¬ 
metry ensures that the operator mixing is the same as 
that in the continuum . Sp ecifically we apply the Rome- 
Southampton method [36| at u = 1.53 GeV, to introduce 
RI/SMOM normalization j2J and then use continuum 
QCD perturbation theory [37| to relate this to the Mini¬ 
mal Subtraction (MS) normalization used for the Wilson 
coefficients Q. 


ANALYSIS AND RESULTS 

The K —> TTTT matrix elements of the operators Qi can 
be determined from the time dependence of the three- 
point functions defined in Eq. 0: 

x(0|J,,(0)K^)(^^|Q,(0)|A)(A|J^(0)|0) + -- - . (3) 

The ellipses represent contributions from the vacuum fi¬ 
nal state or excited kaon or tttt states. For the “split- 
pion” operator is the time closest to tg. 

The normalization factors (0| J.n.^(0)|7r7r) and 
(A|Jii'(0)|0) in Eq. (jS]), and the energies Mk and 
can be determined from the two-point functions: 

(0|4(G)Jx(4)|0) = (4) 

where X = tttt or K. For X = tttt the contribution 
of the vacuum intermediate state to the left-hand side 
must be subtracted. Figure [2] shows the resulting effective 
energy of the kaon and two-pion states in lattice units. 
The kaon mass is obtained from an uncorrelated fit using 
6 < t < 32. For the more challenging / = 0, tttt energy, 
we perform a correlated, single-state fit over the interval 
6 < t < 25, obtaining y^/dof = 1.56(68). A correlated, 
two-state fit using 3 < t < 25 gives consistent results. 
We find Mk = 490.6(2.4) MeV and = 49^1) MeV. 
Using the Liischer quantization condition |38l.l39l| we find 
an / = 0, TTTT phase shift So = 23.8(4.9)(1.2)°, smaller 


than phenomenological expectations [^, l4l|. Here the 
first error is statistical and the second an estimate of the 
O(a^) error. For / = 2 we find = 573.0(2.9) MeV 
and will use S 2 = —11.6(2.5)(1.2)°, a corrected version of 
our continuum result [16|. 



FIG. 2. Effective energies of the kaon (squares) and two- 
pion (circles) states deduced from the corresponding two- 
point functions by equating the results from two time sep¬ 
arations to the function A cosh i5efi(r/2 — t) where T = 64 is 
the temporal lattice size, plotted as a function of the smallest 
of those two separations. (We replace T by T — 8 for the tttt 
case.) These are overlaid by the errorbands corresponding to 
the fitted values of (light blue) and rriK (pink). 


Important for type 3 and 4 diagrams is the quadrat- 
ically divergent quark loop. This contribution is the 
same as that from the operator with a coefficient 
oc (rus — •mi)/a^. Since is the divergence of an ax¬ 
ial current, its matrix element between states with equal 
four momentum will vanish and it will not contribute to 
a physical process such as AT —> tttt. However, for ma¬ 
trix elements between states with unequal energies, this 
term may be 20x larger than the other physical terms. 
Even for an energy conserving amplitude, it will con¬ 
tribute both noise and increased systematic error from 
enhanced, energy nonconserving, excited-state contam¬ 
ination. We determine the size of such an unphysical 
piece from the ratio Vi = {0\Qi{tQ)\K)/{0\d'y^s{tQ)\K) 
and then subtract, time slice by time slice, the operator 
rid'y^sitg) [ 4 ^, dramatically reducing the noise for Qs, 
Qgi Qi and Q g. 

The largest contributions to the real and imaginary 
parts of Ao come from Q2 and Qq, respectively. Fig¬ 
ure [3] shows the three-point functions for these operators 
as a function of the time separation between Qi and 
Because the vacuum state may appear between these op¬ 
erators, the relative size of the statistical noise in the 
vacuum-subtracted matrix element increases rapidly as 
— tg increases. In Fig. [3] we have combined the data 
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(by taking an error-weighted average) from each three- 
point function for fixed ~ tQ ^-nd tg — Ik ^ 6. 
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FIG. 3. The Q 2 and Qe three-point functions, plotted in 
lattice units as functions of with the time dependence 

in Eq. (O removed. The horizontal lines show the central 
value and errors from the fit described below. 


i 

Re(Ao)(GeV) 

Im(Ao)(GeV) 


1 

1.02(0.20)(0.07) X 10”' 

0 


2 

3.63(0.91)(0.28) x 10"^ 

0 


3 

-1.19(1.58)(1.12) X 10““ 

1.54(2.04)(1.45) X 

10““ 

4 

-1.86(0.63)(0.33) x 10“® 

1.82(0.62)(0.32) x 

10““ 

5 

-8.72(2.17)(1.80) X 10““ 

1.57(0.39)(0.32) x 

10““ 

6 

3.33(0.85)(0.22) x 10“® 

-3.57(0.91)(0.24) x 

10““ 

7 

2.40(0.41)(0.00) X 10““ 

8.55(1.45)(0.00) X 

10“14 

8 

-1.33(0.04)(0.00) X 10““ 

-1.71(0.05)(0.00) X 

10““ 

9 

-7.12(1.90)(0.46) X 10““ 

-2.43(0.65)(0.16) x 

10““ 

10 

7.57(2.72)(0.71) x 10““ 

-4.74(1.70)(0.44) x 

10““ 

Tot 

4.66(0.96)(0.27) x 10“' 

-1.90(1.19)(0.32) X 10““ 


TABLE I. Contributions to Aq from the ten continuum, MS 
operators for /r = 1.53 GeV. Two statistical errors 

are shown: one from the lattice matrix element (left) and one 
from the lattice to MS conversion (right). See the Supplemen¬ 
tal Material at [URL to be inserted] for tables of the separate 
matrix elements in the lattice, RI/SMOM and MS schemes 
as well as the renormalization matrices which relate them. 


contributions to Aq as listed in Tab. ID Adding these in¬ 
dividual contributions together gives our final result: 

Re(Ao) = 4.66(1.00)(1.26) x 10”'^ GeV (5) 

Im(Ao) = -1.90(1.23)(1.08) x 10”“ GeV (6) 


where the first error is statistical and the second (dis¬ 
cussed below) is systematic. We can then compute the 
experimental measure of direct GP violation: 



f rimA2 

= Re <-;=- ——— 

I V2e [ReAs 

= 1.38(5.15)(4.59) x 10“^, 


ImAo 

ReAn 


( 7 ) 

( 8 ) 


The matrix elements {(7r7r|(5i| Ar)}i<i<io are obtained 
by fitting the corresponding three-point functions to the 
time dependence in Eq. using tq — tx > 6 and 
— tQ>A. We fit 25 time separations with — tx = 
10, 12, 14, 16 and 18. Figure. [3] is consistent with the ex¬ 
istence of plateaus for > 4 and consistent results 

are obtained when including the t-^Tr — tq = i data, sug¬ 
gesting substatistical, excited-state contamination. We 
estimate the systematic error from excited-state contam¬ 
ination as the 5% difference between the tttt amplitude 
from a correlated, single-state fit to the tttt correlator 
with t > A (our matrix element fitting method) and the 
lowest energy amplitude found in a correlated, two-state 
fit to the same data with f > 3, although the difference 
is again within the now smaller statistical errors. (If we 
omit the accurate, t^^T^—tq = 4 data, our statistical errors 
increase by 40%.) Gombining the data into bins of size 1, 
2, 4 and 8 configurations, shows no bin-size dependence 
of the statistical errors, suggesting that autocorrelations 
can be neglected. We therefore use a bin size of one. 

Finally these lattice matrix elements are combined 
with the renormalization factors, Wilson coefficients and 
Lellouch-Liischer finite-volume correction to obtain their 


obtained using the Im(Ao) and Jq values given above and 
our earlier results for Im(A 2 ) and 62 We use the 

experimental values for Re(Ao), Re(A 2 ) and their ratio w 
(since these are accurately determined from the measured 
K —>■ TTTT decay rates) and the experimental value for e. 

We now briefly describe the systematic error estimates 
given in Tab. In] more complete explanations will ap¬ 
pear in a later paper. We estimate the finite lattice 
spacing error by averaging the differences between the 
three, individual A/ = 3/2, FC — >■ tttt matrix elements 
obtained using the present g aug e action [1^ and our re¬ 
cent continuum-limit results |16| . The errors arising from 
the Wilson coefficients are estimated as the difference of 
our result computed using the leading-order (LO) and 
next-to-leading-order (NLO) formulae for Re(Ao) [sj. A 
similar uncertainty arises when we relate our lattice op¬ 
erators to the MS operators in the continuum expression 
for Hw- This procedure is compromised by our use of 
NLO perturbation theory at /r = 1.53 GeV to relate the 
RT and MS-normalized operators and by our omission of 
dimension-5 and 6 quark-bilinear operators (whose con¬ 
tribution we expect to be small) from the nonperturba- 
tive operator matching. These operator normalization 
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errors are estimated, as in Ref. 1^ , by comparing two dif¬ 
ferent RI/SMOM schemes. Parametric uncertainties are 
found by propagating the standard model input param¬ 
eter errors. Comparing two ansatze for the depen¬ 
dence of (5o suggests a 11% uncertainty in the Lellouch- 
Liischer finite-volume correction. Finally systematic er¬ 
rors are introduced by our mildly unphysical kinematics 
which are estimated from a companion calculation using 
a 10% larger value of the strange quark mass. 


Description 

Error 

Description 

Error 

Finite lattice spacing 

12 % 

Finite volume 

7% 

Wilson coefficients 

12 % 

Excited states 

< 5% 

Parametric errors 

5% 

Operator renormalization 

15% 

Unphysical kinematics < 3% 

Lellouch-Liischer factor 

11 % 


Total (added in quadrature) 27% 


TABLE II. Representative, fractional systematic errors for the 
individual operator contributions to Re(Ao) and Im(Ao). 


CONCLUSION 


We have presented the first calculation of the direct CP 
violation parameter e' with controlled errors. While the 
2 .1(7 difference between our value for Re(eV£) and exper¬ 
iment gives a strong motivation to refine the present cal¬ 
culation, we believe that the absolute size of our statisti¬ 
cal and systematic errors demonstrates that this is now a 
quantity accessible to lattice QCD. Also for the first time, 
we have computed the real part of the decay amplitude 
Aq. The result agrees with the experimental value and 
provides a test of our methods. This result for Re(Ao) is 
consistent with our earlier explanation of the A/ = 1/2 
rule [i^ in which the large ratio of Re(Ao)/Re(A 2 ) re¬ 
sulted from a significant cancellation between the two 
dominant terms contributing to Re(A 2 ), a cancellation 
which does not occur for Re(Ao). We emphasize that 
this calculation can be substantially improved by adding 
more statistics and by studying larger volumes and addi¬ 
tional lattice spacings to better control the large system¬ 
atic errors. Nonperturbative, step-scaling methods can 
relate the lattice operators being used to those defined at 
much smaller lattice spacing where the perturbative Wil¬ 
son coefficients can be more accurately determined. We 
expect that a 10% error relative to the measured value 
of Re(£'/£:) can be achieved within 5 years, motivating 
continued improvement in the experimental result. Sub¬ 
stantially more accurate results will become possible with 
further increases in computer power and the inclusion of 
electromagnetism. 
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i 

Ml:i {GeVf 

{GeVf 

i 

■^SMOM (GeV)^ 

(GeV)3 

1 

-0.247(62) 

-0.147(242) 

1 

-0.0675(1109)(128) 

-0.151(29)(36) 

2 

0.266(72) 

-0.218(54) 

2 

-0.156(27)(30) 

0.169(42)(41) 

3 

-0.064(183) 

0.295(59) 

3 

0.212(52)(40) 

-0.0492(652)(118) 

4 

0.444(189) 

— 

4 

— 

0.271(93)(65) 

5 

-0.601(146) 

-0.601(146) 

5 

-0.193(62)(37) 

-0.191(48)(46) 

6 

-1.188(287) 

-1.188(287) 

6 

-0.366(103)(70) 

-0.379(97)(91) 

7 

1.33(8) 

1.33(8) 

7 

0.225(37)(43) 

0.219(37)(53) 

8 

4.65(14) 

4.65(15) 

8 

1.65(5)(31) 

1.72(6)(41) 

9 

-0.345(97) 

— 

9 

— 

-0.202(54)(49) 

10 

0.176(100) 

— 

10 

— 

0.118(42)(28) 


TABLE S I. The unrenormalized matrix elements in the con¬ 
ventional ten-operator basis (second column) and the seven- 
operator chiral basis (third column). Only statistical errors 
are shown. 


In this section we provide additional technical infor¬ 
mation regarding our results. A complete account of the 
calculation will be published in a separate paper. 


Lattice matrix elements 


Table [ST] shows the K —>■ TTTT matrix elements of the un¬ 
renormalized lattice operators. The second column con¬ 
tains the matrix elements of the ten operators {Qi}i<i<io 
written in the traditional, physical basis defined in Eqs. 
(2)-(5) of Ref. 37|. The third column shows the ma¬ 
trix elements of the seven unrenormalized operators in 
the chiral basis, Q', where j sH, 2, 3, 5, 6, 7, 8}. These 
are defined in Eq. (10) of Ref. 37|. For the convenience 
of the reader, we have converted the values in this ta¬ 
ble into physical units and have applied the Lellouch- 
Liischer finite-volume correction of F = 23.96(61) to ob¬ 
tain infinite-volume matrix elements with the normaliza¬ 
tion conventions specified in Ref. 0. 

Only seven of the ten operators shown in the second 
column of Table 1571 are independent; the remaining three 
are related to those by Eq. (9) of Ref. [37|. As we make 
use of stochastic sources to approximate the high eigen- 
modes of our all-to-all propagators, our matrix elements 
only obey these relations at the percent scale. To con¬ 
vert to the chiral basis we discard the matrix elements 
corresponding to the operators Qi, Qq and Qio. 


Non-perturbative renormalization 

The seven unrenormalized operators Q) in the chiral 
basis are transformed to the operators renormal¬ 
ized in the non-perturbatively defined Rl/SMOM(^, 
scheme, by the 7x7 matrix 

Qf = (SI) 

j 


TABLE S11. The renormalized matrix elements in the 
Rl/SMOM(^, ^) scheme and chiral basis (second column). 
The third column shows th^ matrix elements of the tradi¬ 
tional, physical operators defined in the MS scheme. The 
latter are obtained by applying the 10 x 7 scheme-change ma¬ 
trix factors given in Table ISTVl to the numbers in the second 
column. The left error shown is statistical and the right is 
systematic. For the second column the systematic error is 
a uniform 19% estimate obtained from Table 11 by omitting 
the errors associated with operator renormalization, input pa¬ 
rameters and Wilson coefficients. (This estimate ignores the 
errors associated with our use of an incomplete set of off-shell, 
Rl/SMOM operators, an error which we perturbatively esti¬ 
mate at the 1% level.) The errors in the third column are 
similarly obtained from Table 11 except only the parametric 
and Wilson coefficient errors are omitted, giving a uniform 
24% systematic error estimate. 


for k G {1, 2, 3, 5, 6, 7, 8}. The matrix is deter¬ 

mined by requiring that the operators (after the 
subtraction of dimension-three and dimension-four op¬ 
erators) obey the Rl/SMOM (^, f) normalization condi¬ 
tions specified in Ref. • These are conditions imposed 
on amputated, Landau-gauge, four-fermion Green’s func¬ 
tions which contain the operators and are evaluated 
at off-shell, external momenta that are specified by the 
two momenta pi and p2 ■ In our case pi = ^(0,4,4,0) 
and p 2 = ^(4,4, 0,0) which together satisfy the sym¬ 
metric momentum condition = (pi — P 2 )^ = 

(/i)^, where L = 32 is the spatial extent of the lat¬ 
tice. Here p = 1.53 GeV determines the renormalization 
scale. The resulting matrix elements of these transformed 
Rl/SMOM(^, ^) operators are given in Table ISTIl and the 
NPR matrix in Table ISTlIl 

Finally we use the one-loop perturbative formulae 
given in Eqs. (54), (56), (60) and Table XI of Ref. 371 


to evaluate the matrix Z^^^^ which determines the ten 
traditional MS operators in terms of the seven chiral 
basis operators 




k 




(S2) 


It is the matrix elements of these ten MS operators, mul¬ 
tiplied by the appropriate Wilson coefficients, which are 
given in Table HI 
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0.4582804(281) 

0 

0 

0 

0 

0 

0 


0 

0.2928(300) 

-0.3130(349) 

-0.0123(112) 

0.00579(714) 

0 

0 

0 

0.1368(481) 

0.8350(573) 

-0.0051(209) 

0.0060(139) 

0 

0 


0 

-0.160(134) 

-0.160(143) 

0.4393(486) 

-0.0707(286) 

0 

0 


0 

0.0335(786) 

0.106(100) 

-0.0355(367) 

0.3456(239) 

0 

0 


0 

0 

0 

0 

0 

0.4923552(284) 

-0.0921608(285) 

) 

0 

0 

0 

0 

0 

-0.0593122(599) 

0.372432(125) 


TABLE Sill. The 7x7, lattice-to-RI/SMOM(^, conversion matrix Only statistical errors are shown. 


/ 0.19744 

1.07389 

0.13830 

0 

0 

0 

0 

\ 

0.19744 

-0.12667 

0.75559 

0.00208 

-0.00625 

0 

0 


0 

2.96832 

1.92607 

0.00417 

-0.01250 

0 

0 


0 

1.74719 

2.49297 

0.01554 

-0.04636 

0 

0 


0 

0 

0 

1.00121 

-0.00364 

0 

0 


0 

-0.04687 

-0.10936 

-0.00164 

0.99519 

0 

0 


0 

0 

0 

0 

0 

1.00121 

-0.00364 


0 

0 

0 

0 

0 

-0.01726 

1.04206 


0.29616 

0.12667 

-0.75559 

-0.00208 

0.00625 

0 

0 


V 0.29616 

-1.07389 

-0.13830 

0 

0 

0 

0 

/ 


TABLE SIV. The 10 x 7 RI/SMOM(^, —>■ MS conversion matrix, computed to one-loop in QCD perturbation 

theory at the scale p = 1.531 GeV. 


1 

Zi 

Vi 

1 

-0.3734 

0 

2 

1.189 

0 

3 

0.001294 

0.02720 

4 

-0.003691 

-0.05828 

5 

0.002505 

0.007133 

6 

-0.004718 

-0.08191 

7 

6.099x10"® 

-0.0003401 

8 

-4.414x10"® 

0.0008684 

9 

3.561x10"® 

-0.01046 

10 

3.010x10"® 

0.003496 


TABLE SV. The Wilson coefficients in the MS scheme at 
At = 1.531 GeV. 

Wilson coefficients 

The real and imaginary parts of the amplitude Aq = 
{{'mT)i=Q\Hw\K) are obtained by multiplying the MS- 


renormalized matrix elements by the appropriate Wilson 
coefficients and CKM matrix elements as given in Eq. O- 


The Wilson coefficients were computed at the renor¬ 
malization scale of /r = 1.531 GeV from the equations 
given in Ref. @. For the required standard model param¬ 
eters we use the values listed in Table IS Vli from which 
we obtain the Wilson coefficients shown in Table ISVl We 
use the two-loop beta function to obtain the three-flavor 
value of as = 0.353388 at /r = 1.531 GeV. The GKM 
matrix elements V^s and Vudi along with the remaining 
inputs necessary to compute Re(e'/e) from Eq. ([7]) are 
also listed in Table [SVD 
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Input 

Value 


1.531 GeV 

mw 

80.385 GeV 

mt 

160.0 GeV 

rric 

1.275 GeV 

mb 

4.18 GeV 

^EM 

127.940 

sw?{6w) 

0.23126 

as{mz) 

0.1185 

A 4 

0.331416 GeV 

T 

0.001543 - 0.000635i 

V^s 

0.2253 

Vud 

0.97425 

UJ 

0.04454 

kl 

0.002228 

<t>e 

0.75957 rad 

R'6(-Ao)expt 

3.3201x10"’^ GeV 

R'G(-^2 )expt 

1.479x10"® GeV 


TABLE S VI. Inputs for the computation of the Wilson coefficients are given above the break, and the remaining inputs that 
were used for cormuting Re(£'/£) are given below. The parameters given in this table were obtained from the PDG Review of 
Particle Physics [J. Here (j>e is the phase of the indirect CP-violation parameter e. The last two entries are the experimental 
values of Re(Ao) and Re(A 2 ), which we combine with our lattice values for the imaginary components to obtain Re(£'/e). 
Other than the systematic errors listed in Table HIl we neglect the errors on these inputs as their contributions to our final error 
are small in relation to our statistical errors. 







